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1. INTRODUCTION 

1.1. Mass Spectrometers for Chemical Process Control: 

Process control in the chemical industries is effected 
through various methods. But most of them are confined to 
the control of process parameters like flow rates, levels, 
temperature, pressure etc. While ij is necessary to keep 
these process variables controlled, it is imperative that 
methods of control through product analysis should be evolved 
as the final product quality is of greatest concern to us. 
This implies that procedures of analysis of final products 
should be evolved and the results of this analysis should be 
transformable into controlling signals. Some methods of 
analysis have come into vogue in the past such as refract o- 
metry, colorimetry, absorption spectrometry, polorography, 
gas chromotography etc. All these methods are based on the 
physical properties of the sample such as mechanical proper- 
ties, electrical properties, thermal properties etc. The 
use of these methods remain limited due to -this dependence. 

A recent entry into this field has been made by mass 
spectrometers which do not have this disadvantage of the 
other analytical instruments. Mass spectrometers determine 
directly the mass number of the sample species. Currently a 
combination -of gas chromotographs with mass spectrometers is 
widely used for chemical analysis. Of the great variety of 
mass spectrometers available, Quadrupole mass spectrometers 
offer some advantage over others because of their character- 
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i sties. Since these are also less expensive Quadrupole mass 
spectrometers have the potential of in situ process control 
in chemical industries. 

1.2 Basic operation of Quadrupole Mass Spectrometer i 

Unlike most other mass spectrometers, the quadrupole 
mass spectrometer does not require a magnetic field. A 
typical quadrupole mass spectrometer is shown in Pig. 1, 

The principal elements constituting the mass spectrometer 
are the ion source, the quadrupole mass analyser and the ion 
collector (or detector). The whole system operates under 
high vacuum. 

The gas to be analyzed is ionised by bombardment by an 
electron beam which is usually from a hot filament. As a 
result of this electron bombardment a small percentage of the 
gas atoms or molecules gets ionized. These ions are accelera- 
ted and through an entrance aperture enter the radio- frequency 
electric field in the mass analyser region. The field 
conditions (to. be discussed later) cause a filtering action on 
the ion beam. This filtering results from the mass analyser 
permitting only those ions of a specific range of mass to 
charge ratio to completely pass the field length. These ions 
have bounded trajectories and they pass through to the exit 
aperture. All other ions have unbounded trajectories and are 
collected on the rods constituting the quadrupole. 

The ions that travel through the analyser are ulti- 
mately collected- on the ion collector. The output current 
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from the collector is proportional to the amount of the 
particular species in the given sample, assuming idealized 
conditions in the ion source. The charge to .mass ratio 
satisfying the criterion for a hounded trajectory and' passage 
through the field length is linearly related to the potentials 
on the quadrupole rods. Thus calibration with any known pure 
gas sample is sufficient to identify any other component-. 

1.3 Advantages of the Quadrupole Mass Spectrometer: 

Since the quadrupole mass spectrometer is non-magnetic , 
the instrument is light in weight and problems of magnetic 
screening where stray magnetic fields might interfere with 
the system under observation are avoided. Since the principle 
of mass separation in the quadrupole mass analyser is un- 
affected by the collisions between the particles, the instru- 
ment can be used at relatively high pressure. The mass 
analyser’s relative insensitivity to the initial velocity of 
ions along the analyser axis enables its use where there is 
a spread in the initial ion energy. The linear relation 
between the charge to mass ratio and the quadrupole electric 
potential, makes it simpler to use. The mass spectrometer is 
simple in' its configuration and compact in design. These 
features of this mass spectrometer permit its wide use in 
various fields. 

1,4. Objectives of this Work: 

s^nln this work the per mis sable entry conditions of ions 
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have been computed. These data aid in designing the ion 
source and entry aperture of the mass analyser. A medium 
resolution mass spectrometer has been fabricated and has been 
used as a residual gas analyser. ^ 

*** 
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2. THEORY OF OPERATION OF QUADRUPOLE MASS 

SPECTROMETER 

2,1, Electrode Structure arid Quadrupole Field: 

The quadrupole electric field is such that it is 
linearly related to the distance of the ion from the axis. 

This type of field has the nature of focusing the ions at the 
axis. This field is produced by a potential of the form 

0 = { t Q (ax 2 + by 2 + cz 2 ) ... (1) 

where a,b,c are constants characteristic of the particular 
field. If there is no space charge within the electrode 
structure, the potential must satisfy Laplace’s equation 
2 0 = 0. This yields 

a + b + c = 0 ... (2) 

If the field is two-dimensional and taking z to be the axis 
of symmetry, we get 

c = Oj a = -b = l/ r £ (say) ... (3) 

d = d (£%C) ... (4) 

r 0 

This potential can be produced by a configuration as in Fig. 2. 

of four hyperbolic surfaces located so as to correspond to 

equipotential lines of the quadrupole field. The voltage 

0 q = (U + V cos wt) is applied to the electrodes with the 

centre of the field chosen' as the zero of the potential. The 

potential at any point is then 

2 2 

0 = (U + V cos vrt) ... (5) 

r o 

The x direction is that through the centres of the pair of 
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rods for which the dc component is positive with respect to 
the second pair, 

2.2 Equations of Ion Motion: 

The electric fields obtained from the quadrupole po- 
tential (Eq. 5 ) are 


E y ■ 


= 


d 0 

3x 

M 

O0 

2 > z 


-(U + Tcoswt)~ 

*o 

+ (U + 7 cos wt) 

r o 


0 


( 6 ) 


The equations of motion of a singly charged ion of mass m are 
2 

d 0 f TT Oy 

+ 5 (U + V cos wt) ~ =0 

r o 

2 

7-f - | (U + 7 COS wt) = 0 ... ( 7 ) 

d 2 z/dt 2 10 = 0 

The motion in each coordinate direction is thus independent 

of the motion in the other directions. Introducing the substi- 
tutions 

wt = 2 £ 

a = 8eU/mr 2 co 2 
q = 4 eV/mr 2 60 2 
(Bq. 7 ) can be reduced to 

2 o 

d x/dfef + (a 4- 2c[ cos 2^ 

^ y/ds£, 2 - (a + 2q cos 2^ 

d 2 z/d^ 2 
— ^ 

These differential equations are of the type known as Mathieu's 


( 8 ) 


) x = 0 

) y = o (9) 

= 0 
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equations. 3? lie properties of the Mathieu equations are well 
known(l). We are primarily interested in the solutions that 
result in stable or bounded ion motion. 

2.3 Solutions to the Equations of Motion: 

Solutions to the Mathieu equation can be expressed in 
the form 

L 4-'0 t 

+ *c"e-^C zs e- Zi -°klO) 

where u is x or y. The solutions are of two types, solutions 

denoted "stable" where u remains finite as and "unstable-* 

where u increases without limit as^a^. The constant 

determines the type of solution.,/! depends upon the values of 

a and q, and is independent of the initial conditions. If /*- 

is real and non-zero, instability will arise from either the 
-M-k 

e 7 or e factor. There are three other possibilities: 

1. j* - iffy is purely imaginary and . jjb is not a whole number. 

The solutions are stable. 

2. ja is a complex number. The solutions are unstable (assum- 
ing that u Q and u^ are not both zero) 

3. /A. = in is purely imaginary and n is an integer. The 
solutions are periodic but unstable. 

Since /a, depends on a and q, the regions of stable 
operation can be plotted in a-q space. Such a figure is known 
as the stability digram. The stable areas with large a, q 
values correspond to ion motion in which the amplitudes of 



10 


oscillation are too- large to be useful in practical devices. 
The stability area near the origin is of practical signifi- 
cance. The stability limits are (p> = 0 and (ft = 1. 

2.4 Stability Diagram: 

Dor an ion to traverse the mass filter, it must have 
stable trajectories in both x and y directions. We can super- 
impose the stability diagrams for x and y directions to deter- 
mine the stability region for the mass filter as shown in 
Dig. 3. A line on the stability diagram with slope a/q called 
mass scan line is defined if the ratio of dc to rf voltage is 
fixed since a/q e'quals 2U/V. The apex of the stability tri- 
angle lies at q-j . . t = 0.70600 and = 0. 23699. Consider 

a mass scan line as shown in Dig. 3. Dor any chosen values of 
U and hence V, ions of differing e/m will be spread out along 
the line according to their (a,q) values, with ions of larger 
e/m farthest from the origin. This mass spread can be reduced 
by increasing the ratio a/q i.e. U/V. This spread can be made 
as narrow as desired to permit only one mass to pass through 
the field. All the other ions will execute unstable oscilla- 
tions, strike the electrodes and will be removed. This is the 
fundamental filtering action of the quadrupole field. 

The mass scanning can be done by varying any of the 
field conditions^, U or V since practically r Q cannot be 
varied. Maintaining tuned circuits in the amplifier stages 
over a frequency range is the most serious drawback of the 
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variable frequency method. Since a and q vary as — ^ the 
frequency would also have to change as the square root of the 
control function if the mass peaks are to appear in a linear 
fashion. Hence varying CD is generally not recommended. 

Keeping either one of them constant, U or V can be changed to 
obtain the mass spectrum. But this results in shifting the 
scan line during the scanning thus introducing a continuous 
variation in resolution. Hence this method is also not re- 
commended. By varying the magnitude of U or V keeping constant 
the ratio U/V and CD , each ionic species may be brought in turn 
into the area for stability in both x and y directions and can 
be transmitted through the device. Ions with e/m greater than 
the stable species have unstable trajectories in the x direction 
and ions with e/m lower than stable ions are unstable in y 
direction. 

2.5 Ion Trajectories in the Quadrupole field: 

Typical ion trajectories in the xz and yz planes are 
illustrated in fig. 4 and fig. 5 respectively. The stability 
criterion of ions in the field does not depend on the initial 
conditions of ions at the entrance aperture such as the initial 
velocity (radial or axial), initial position and entrance phase 
angle. However this independence is somewhat hypothetical 
since the stability criterion considers only the conditions that 
the ions be either bounded or unbounded. The initial conditions 
have an effect on the practical stability since the field 
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i.e. field diameter .and field length, are finite. 

If the radial velocity of ions at entry exceeds a 
limiting value, then even though the ions may have hounded 
oscillations, they would strike the electrode and he removed. 
Hence there would he a reduction in the expected value of 
transmission of ions. If the radial velocity of the ions 
at entry is helow a limiting value, then even though the ion 
oscillations may he unstable it would pass through the finite 
field length without striking the electrodes. Hence there 
would he a reduction in the expected value of resolution. 

Similar limits are imposed hy the initial position of ions also. 
Hence of practical importance here, is the placing of limits 
on the initial displacements x q and y , and the initial radial 
velocities and y Q . The axial velocity at entry is important 
to the extent that . it must he of low enough magnitude to 
ensure that the ion stays in the electric field long enough to 
separate the ions of adjacent mass numbers for a given resolu- 
tion. An increase in resolution i.e. an increase in the U/V 
ratio has the effect of increasing sharply the amplitude of ion 
oscillation. Thus for given initial conditions an increase in 
resolution i.s accompanied hy reduced ion transmission. So, 
in practice, a balance between resolution and transmission 
should be struck. 
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3. EFFECT OF INITIAL CONDITIONS ON ION TRANSMISSION 

3*1 Maximum Amplitude of Ion Oscillation: 

It was shown earlier that in a practical device, not 
only must a confined ion have a stable trajectory, but also the 
maximum amplitude of this trajectory must be smaller than the 
field radius. This necessitates computation of maximum ampli- 
tudes of ion motion as a function of the initial displacements 
x Q and y Q , initial radial velocities x^ and y^ and initial 
phase cot Q for (a,q) values chosen in the region of interest 
in the stability diagram. The equations of motion were numeri- 
cally solved using appropriate initial conditions x 0 ,y ,x^,y^ 
and u) t_ to determine the maximum amplitudes of ion trajectories. 
Fourth order Runge-Kutta method was- used for this purpose. 
Computations were made for various operating points along the 
mass scan lines corresponding to resolutions in the range 70— 
1000. The maximum amplitude of ions was determined over 160 
cycles of ion motion. 

3.2 Results of Computation: 

The ratio of maximum amplitudes of ions to initial 
x y 

position — and ~ are plotted against * q f values for reso- 
x o y o 

lutions 100,200,500 and 1000 in Fig. 6. The f a' values at 
which the computations are done are those corresponding to the 
given resolutions. The initial position of ions is x o =y Q =0.01 
inch (the amplitudes can, however, be arbitrarily scaled) and 
the radial velocity at entry is x^ = y^ = 0.005 inch/rad. The 
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The initial phase ,k*t is zero. 

The effect of change in initial radial velocity of 
ions for a resolution of 200 is shown in Fig. 7, The initial 
displacement of ions is 0.005 inch and the initial phase at 
entry is zero. The initial .radial velocities are 0.005, 0.0075, 
0,0025 inch/rad. respectively. Similar results for a reso- 
lution of 500 are shown in Fig. 8. 

Fig. 9 shows the effect of change of initial position 
at a resolution of 70. The radial velocity of ions at entry 
is 0.01 and phase at entry is zero. 

The effect of initial phase on the maximum amplitude 
of ion motion are shown in Fig. 10, for operating points along 
different iso-^> lines. The initial position of ion at entry 
is 0.01 inch and the radial velocity at entry is 0.005. 

Using the above plots the limiting values of ' q* where 

\ 

x and y amplitude equals field radius is determined for given 
initial conditions. Within this range of 'q' all stable ions 
are transmitted resulting in 100$ ion transmission. These 
ranges for 100$ ion transmission for given initial conditions 
of ions are represented by contours in Fig. 11. The field 
radius is taken to be 0.350 inch to arrive at these 100$ ion 
transmission contours. 

Max imum amplitudes of ion motion computed for parameter© 

x , x\ cot and A given by Paul (2) agree within 1$ to Paul’s 
0*0*0 V 

computations (2), 
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3.3 Conclusions from Present Computations: 

Considering the effects of initial position of ions 

and radial velocity at entry independently (i.e. either of 

them is zero while the other is non-zero) Paul (2) arrived 

at limiting values for x and x as 

oo 

2y 0 , 2x 0 ~ r o ( hM/M) y 2 ... < 1 1) 

V <C °‘ 16 r o W ( 2^ M/M) ^2 ... (12) 

for 100$ ion transmission, The present computations with both 
initial conditions non-zero simultaneously yields the bounds 
given in Pig. 1 1 for 100$ ion transmission which are higher 
than the bounds given by Eq.11. and Eq. 12. This is presu- 
mably due to the fact that these equations do not account for 
interaction effects when both initial conditions are non-zero. 
Brubaker (3) has experimentally observed that off-axis entry 
of ions (i.e. non-zeyo initial displacement) with non- zero 
initial radial velocity results in higher ion transmission 
than ion entry at axis (i.e. zero initial displacement). This 
implies that interaction effects (i.e, simultaneous non- zero 
initial conditions) result in higher bounds for 100$ ion 
transmission than individual effects (i.e., either of the 
initial conditions is zero). These computed values are used 
in determining the limits of initial conditions for 100$ ion 
transmission at given' resolution. 
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4. QUADRUPOLE MSS SPECTROMETER DESIGN 
4.1, Mechanical Design: 

'The ion source is of electron impact type. A schematic 
diagram is shown in Dig. 12. The electrons are emitted by a 5 
mil diameter heated Tungsten filament. An accelerating 
potential applied between the filament and the ionisation 
chamber takes the electrons into the ionisation chamber. The 
sample to be analysed, from a sampling device, enters the 
ionisation chamber via a passage in the repeller. Electron 
impact ionisation takes place within the ionisation chamber 
and the electrons resulting from ionisation are removed at the 
anode which is on the other side of the ionisation chamber 
opposite to the filament. The positive ions produced in the 
ionisation chamber are repelled away by the repeller and enter 
into the mass analyser section through an entrance aperture 
of 0.020 inch diameter. 

The mass analyser consists of four cylindrical rods 
in place of hyperbolic rods. They are arranged such that 
r rod //r o = (2). This configuration results in best 

approximation of the hyperbolic field using cylindrical rods. 

The quadrupole rods are 6" long and 3/8 n in diameter 
and are machined accurate to 0.001 inch. The end supporting 
pieces on which these rods are mounted are obtained bycutting 
a machined piece into two parts. The mounting holes are bored 
and the alignment channels are cut before it is separated into 
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two end pieces. The mounting of quadrupole rods is shovm in 
Pig. 15. The alignment of the rods is accurate to 0.001 inch. 
The rods in the field region are insulated from the rod pieces 
at support by special ceramic insulators. These ceramic insu- 
lators are metal coated at their, ends and are silver soldered 
to the quadrupole rods. By this way, the tolerance of quadru- 
pole rod alignment is made independent of insulator non- 
uniformities. 

The ion collector used at the end of the quadrupole 
field is an EGA electron multiplier. It is a 14 stage Gu-Be 
dynode multiplier with a gain of 10^. 

The whole assembly of ion source, analyser and electron 
multiplier ion collector is flange mounted. It is compact, li 
light in weight and can be easily mounted on any system. 

4.2. Electronic Design; 

The electronic circuit to operate the quadrupole mass 
spectrometer is given in Pig. 14 and 15. The 1 MHz rf voltage 
is produced with a crystal oscillator and power amplifications 
of this rf voltage is accomplished in the tube circuit shown 
in Pig. 14. The dc voltage is obtained by rectifying a part 
of the rf voltage taken from an inductive voltage divider 
shown in Pig. 15. By this the constancy of the ratio between 
rf and dc voltages is ensured. The inductive voltage divider 
is grounded at centre thus symmetrising the field. Two 
potentiometers on either sides of this ground are used to vary 
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the ratio U/V. 'The frequency is staple to 10"”^. But consi- 
derable difficulty is encountered in symmetrising the field. 
The mass scanning is accomplished by varying the screen 
voltage to the oscillator amplifier tube. 

Big. 1 6 gives an alternative transistorised circuit 
which was also used. This transistorised rf amplifier is 
used in conjunction with an oscillator (output voltage of 
0-5 V) and the AG/DC matching network shown in Big. 15. When 
this is used as the quadrupole power supply, the mass scanning 
is accomplished by varying the oscillator output voltage. 
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5. EXPERIMENTS. RESULTS AND CONCLUSIONS 

The quadrupole mass spectrometer was tested with, gas 

mixture containing N £ , 0 2 , Ar, HO, C0 2 * The operational 

pressure was about 10 mm Hg measured by an ionisation gauge. 

The predominant ions under these experimental conditions are 

^2 * ^ + > -^ r+ » N0 + , 0+ 0 + , 00+ The unanalysed ion current 

-“11 f 

is of the order of 10 ampere using 10 as the gain of the 
electron multiplier detector. The variation of the ion current 
is linear with pressure within experimental errors. Fig. 17 
indicates some typical mass spectrums obtained with a sample 
of mixture of gases. The observed resolution of the mss 
spectrometer is significantly lower than the design value of 
500. 

For process control applications the ion current can 
be amplified and this amplified signal can be used as the 
controlling signal. The present flange mounted design of the 
mass spectrometer is ideal for its use as in situ controller. 

Since the observed resolution of the spectrometer is 
significantly lower than the design value, further work is 
necessary to improve the resolution. As the present design 
of quadrupole mounting is close to the requirements, any 
increase in resolution can be achieved only with improvements 
in electronic design. Specifically, the voltage balance 
between either pair of the quadrupole rods should be improved 
and the constancy of the U/V ratio should be maintained. These 
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two. improvements should significantly increase the resolution. 
Also, developmental work for the fabrication of electron- 
multiplier particle detectors is to be taken up so that the 
RCA multiplier being used at present can be replaced by our 
own indigeneous design. 
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THIS PROGRAM COMPUTES THE MAXIMUM AMPLITUDES OF AH I ON OSCILLATIONS 
IN THE -X- AND » Y« DIRECTIONS OF THE QUADRIPOLE FIELD FOR VARIOUS 
ENTRY CONDITIONS OF IONS SUCH AS POSITRON* RADIAL VELOCITY AND PH ASF 
AT DIFFERENT MASS SPECTROMETER OPERATING POINTS. 




s mmm® m» ik.-* mm mm mmmm w* mmwmmmmmmm 


I8FTC MAIN 

DIMENSION AK( 12) » AL?12 ) ,F I ( 20 ) , Z TEMP { 500 ) ,FTEMP ( 500) , XT"MP ( 5 00) 
DIMEN5I0NA8I4) 

INDEX=0 

1=1 

H=0. 05 


ABU )=0. 23627 AND ABU ) = 0.23667 CORRESPOND TO THE A( 0=0.706) 
VALUFS FOR GIVEN MASS SPECTROMETER RESOLUTIONS 


ABU ) =0.23627 
AB ( 2 ) =0.23663 
D0107 1 J=l* 2 
AA=AB { I J) 

PR? NT 111. AA 

111 FORMAT (// ?X* ♦AA=*F’* 5.8//) 

R ATIO = AA/O.T 06 
0=0.696 

105 0=0+0.003 

A = C*R ATI 0 
99 FI ( I ) =0.0 
K = 0 

MAX =10 
C 
C 

C ZET A = 0 . 0 CORRESPONDS TQ PHASE AT ENTRY 
C XI N= 0.005 CORRESPONDS TO POSITION AT ENTRY 
C FIN 1=0.0 05 RADIAL VELOCITY AT ENTRY 

C 

c 

c 

59 ZETA=0 . 0 
X I N=0 . 006 
FINI =0.005 
X I N I =X I N 
X=X IN 
FIN=FINI 
F = F IN 

PRINT lOiXINf FIN*A»0 
10 FORMAT (4F16.8 ) 

7 N = 0 
? M=N+1 

ANC*LE=2 »0*ZET A+FK U 
PI=4.0*ATANU .0) 



20 TF5T=ANGLE-2.0*PI 
I F (TEST } 77. 17 , 3.9 

19 AMGLE=TEST 
GO TO 20 

IT PIHLF=ANGL E-PI/2. 0 
IF(PIHLF)21, 21,22 

21 Y=COS( ANGLE) 

GO TO 27 

2? P I0NE= ANGLE- PI 
I F ( PI ONE )23,23,24 

23 8=P I- ANGLE 
Y=— COS ( B ) 

GO TO 27 

24 P IT HR = ANGLE" (3. G*P T ) /2 .0 
IF ( PITHR )25»25»26 

25 Z= ANGLE-PI 
Y=-COS ( Z ) 

GO TO 27 

26 Z=2.0*PI"ANGL£ 

Y=COS ( Z ) 

27 G=-( A+2.0*Q«Y )*X 

I F ( INDEX™! 1 10 0, 101 , , 01, 

101 G*-G 
100 AK (N)=H*F 
AL f N ) =H*G 
X-XIN 
F=FIN 

IF ( N—3 3 5 , 5 

4 ZETA=ZETA+H/2.0 
X=X+AK ( N ) / 2 .0 

F = F+AL ( N ) / 2 .0 

IF{N~2)3,6,6 

6 ZFTA=ZFTA“H/2 .0 
GO TO 3 

5 IF(N“A)7»8 ,8 

7 ZFT ft*ZET 4+H/2 .0 
X = X + A K { M 3 

F = F+AUN) 

GO TO 3 

8 AKSUM* ( AK( X 3 + 2. 0*AK f 2 3 + 2 . Q*AK( 3 3 +AK (433/6.0 

9 ALSUM= ( AL ( 1 )+?.0*ALf 2) +“ . 0*AL { 3 3 +AL< 4) 3 /6.0 
X=X+AKSUM 

F*F+ALSUM 
I F ( X IN 337, 38, 38 

37 IF (X >41*64 *64 

38 IF (X)64,*'l ,41 
41 IF ( M AX 3 46 ,45,45 

45 IF (X-XIN 548,48,64 

46 T F ( XI N) 64, 4-7, 47 
48 MAX*-I 0 

GO TO 51 

47 MAX *1.0 
51 K = K +1 

XTEMP ( K) =X IN 
ZTFMP(K3=ZETA«H 
FTFMP ( K)*F IN 
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64 X I N=X 

FIN = F 
X = X I N 
F = FIN 

I F ( ZET A"80 - 0) 2,11,1’ 


ZTFMP IS THE PHASE OF OSCILLATION WHEN THE MAXIMUM AMPITUDE IS XTFMP 
FTFMP IS THE RADIAL VELOCITY WHEN THE MAXIMUM AMP LITUDF IS XTFMP 


ll PRINT82, {ZTEMP(L),F‘ r EMP(L),XTEMP{L),L=l,K ) 

82 FORMAT (9?'! 2.5 ) 

0049KM = ?.f 500 ■ 

ZTEMP ( KM)=0.0 
FT CMP ( KM 1 = 0*0 
49 XTFMP ( KM )=0-0 

TF( INDE5M) 10 2,?9,T« 

302 INDEX=INOEX+l 
GO TO 99 
T9 INDFX=0.0 

IFIC-0.T08 >105,107, ’07 
1,07 CONTINUE 

900 CONTINUE 

901 FORMAT (2F15. 5 ) 

STOP 

7 : ' END 

£FNTPY 



